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Abstract 

In this talk, we discuss one of the dissipative processes which likely 
take place in the Early Universe. We assume that the matter filling the 
isotropic and homogeneous background is to be described by a relativistic 
viscous fluid characterized by an ultra-relativistic equation of state and 
finite bulk viscosity deduced from recent lattice QCD calculations and 
heavy-ion collisions experiments. We concentrate our treatment to bulk 
viscosity as one of the essential dissipative processes in the rapidly ex- 
panding Early Universe and deduce the dependence of the scale factor 
and Hubble parameter on the comoving time t. We find that both scale 
factor and Hubble parameter are finite at t = 0, revering to absence of 
singularity. We also find that their evolution apparently differs from the 
one resulting in when assuming that the background matter is an ideal 
and non-viscous fluid. 
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1 Introduction 

Dissipative processes are supposed to play a very important role in the Early 
Universe. The most essential ones are the bulk and shear viscosity. First at- 
tempts to create a theory of relativistic fluids, which likely simulate the matter 
under the extreme conditions of the Early Universe, have been done by Eckart 
[1] and Landau and Lifshitz [2]. Regardless the choice of the equations of state, 
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all equilibrium states in these theories (first order) are unstable and in addi- 
tion signals may propagate through the fluid at velocities exceeding the speed 
of light. Therefore these theories are applicable only to phenomena which are 
quasi-stationary, i.e. slowly varying on space and time scales. 

A relativistic "second order" theory has been suggested by Israel [3] and 
developed by Israel and Stewart T] , Hiscock and Lindblom [5] and Hiscock and 
Salmonson [6] into what is called "transient" irreversible thermodynamics. In 
this model deviations from equilibrium (bulk stress, heat flow and shear stress) 
are treated as independent dynamical variables, leading to a total of unknown 
fourteen dynamical fluid variables [7J [5] . 

Causal bulk viscous thermodynamics has been used to describe the dynamics 
and evolution of the Early Universe. But due to the complicated character of 
the evolution equations, very few solutions of gravitational field equations are 
known in the framework of full causal theory. For a homogeneous Universe filled 
with a full causal viscous fluid obeying the relation £ ~ p 1 / 2 [9], with p being 
the energy density, exact general solutions of the field equations have been ob- 
tained in [TU]-[I1]- Because of serious technical reasons, most investigations for 
the dissipative causal cosmologies have assumed Friedmann-Robertson- Walker 
(FRW) symmetry or - in few cases - small perturbations around it [16j . 

Recent RHIC results give a strong indication that in the heavy-ion collisions 
experiments, a hot dense matter can be formed [T7J, which likely agree with 
the "new state of matter" as predicted in the Lattice QCD simulations [IB] , 
According to recent lattice QCD simulations pjjl [20], bulk viscosity £ is not 
negligible near the QCD critical temperature T c . It has been shown that the 
bulk and shear viscosity at high temperature T and weak coupling a s , £ ~ 
alT 3 / \nct~ 1 and rj ~ T 3 /(allnaj 1 ) [21]. Such a behavior obviously reflects 
the fact that near T c QCD is far from being conformal. But at high T, QCD 
approaches conformal invariance, which can be indicated by low trace anomaly 
(e — 3p)/T 4 [22], where e and p are energy and pressure density, respectively. 
In the quenched lattice QCD, the ratio (,/s seems to diverge near T c [25] . 

To avoid the mathematical difficulties accompanied with the Abel second 
type non-homogeneous and non-linear differential equations, one used to model 
the cosmological fluid as an ideal (non- viscous) fluid. No doubt that the viscous 
treatment of the cosmological background should have many essential conse- 
quences |24j . The thermodynamical ones, for instance, can profoundly modify 
the dynamics and configurations of the whole cosmological background [8] . The 
reason is obvious. The bulk viscosity likely arises due to increasing entropy and 
can be expressed as a function of the overall energy density p [S]. 

In this talk, we assume that the background corresponding to FRW model 
is filled with ultra-relativistic viscous matter, whose bulk viscosity and equation 
of state have been deduced from recent heavy-ion collisions experiments and 
lattice QCD simulations. The relativistic fluid will be described according to 
the Israel- Stewart (IS) theory. 
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2 Model 



We take into consideration natural units, i.e., k — c — 1 and h — 2tt and assume 
that background geometry of early Universe is filled with bulk viscous cosmic 
fluid, which can be described by a spatially flat FRW metric. The line element 
is 

ds 2 = dt 2 - a 2 (t) [dr 2 + r 2 (d9 2 + sin 2 6d4> 2 )] . (1) 

At vanishing cosmological constant A, Einstein gravitational field equations in 
a static and flat Universe (k = 0) read 

Rik- -9i k R = ^GT ik . (2) 

Inclusion of bulk viscous effects can be generalized through an effective pressure 
IT, which is formally included in the effective thermodynamic pressure p e f / [7] . 
Then in the co-moving frame the energy momentum tensor has the components 
T ° = p, 2f = T 2 2 = ^3 3 = —Peff- For the line element a given by Eq. (p}, we get 
two independent solutions 

-=-^G (3p o// + p), (4) 

where one dot denotes differentiation with respect to the comoving time t and p 
is the energy density. Assuming that the Universe is a closed system, the total 
energy density of cosmic matter is conserved, i.e., = 0. 

p + 3H( Peff +p) = 0, (5) 

Here we introduced the Hubble parameter H = a/a. In the presence of bulk 
viscous stress n, the effective pressure becomes p e ff = p + H, where p is the 
thermodynamic pressure of the cosmic fluid. Then Eq. ([5]) can be written as 

p + 3H (p + p) = -3UH. (6) 

According to the causal theory of relativistic fluid; Israel-Stewart theory [31(3], 
the evolution equation of the bulk viscous pressure reads [7] 

rtl + U = -^H-^Tu(^H+^-^-^y (7) 

where T is the temperature, £ denotes the bulk viscosity coefficient and r stands 
for the relaxation time. 

Equations of state for p and T can help to have a closed system from Eq. ([3]) , 
([6]) and ([7]). r and £ are determined according to phenomenological approaches. 
For instance, bulk viscosity of QGP at high temperature can be determined 
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approximately from the equilibrium thermodynamical equation of state taken 
from recent lattice QCD calculations [T9l 125) . 

P = ^P, T = (3p r , t = ap+—T2, (8) 

with u = (7 - 1), uj ^ 0.5 - 1.5 GeV, (3 = 0.718, 7 ~ 1.318, r = 0.213 and 

1 97 2 - 24 7 + 16 



a 



9uj 



7-1 



(9) 



From Eq. and ® we obtain an equation describing the cosmological evolu- 
tion of Hubble parameter H 

H + l[l + (l- r)i\HH + -H-(l+r)^- + %- 2)H 3 + —H 2 = (10) 
2 a H 4 a 



With the transformation u = H, Eq. [TO] reads 



u~(l+r)% + + (1 - r ) 1 ]H+ I) u+% -2)^ + 11^ = (11) 



a 

H 



2a 



which is a second kind Abel type first order ordinary differential equation [26] , 
We take into consideration some numerical approximations 27, 28J in order 
to reduce Eq. QT] to an Abel canonical form. We start with the general form of 
the second kind Abel differential equation 



[y + g(x)]^ = f 2 (x)y 2 + h(x)y + f (x), 



(12) 



where g(x) 7^ 0. Last equation can be reduced to the Abel canonical form [26] 
by using the following transformations 



[y + g(x)]E 



e = 

Then the canonical form reads 



exp 



f2(x)dx 



where, 



-— = F 1 (x)u + F (x) 
ax 

F^x) = (f 1 -2f 2 g + gQE 
Fq(x) = (/o - fig + f 2 g 2 )E 2 



(13) 



Introducing new independent variable z = J Fi(x)dx in Eq. (|13| results in 

uiuj' z -u) = g{z) (14) 
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where g(z) is defined parametrically as [271 128] 



g(z) = ^^Az + B (15) 

where A and B are constants. For simplicity, we approximate g{z) to A z. Then 
Eq. HU turns to be solved in H (t) 

H(t) = (16) 

where denominator is dimcnsionlcss, 

-3[1 + (1 - r) 7 ] a- 1 

2(1^ 7) ' ~> £ = exp(-Bat) 

The dimcnsionlcss a is taken as a free parameter. Obviously, it's physical value 
in the expanding Universe is likely negative. Finally, the scale factor reads, 



a{t)=aa {'£~A) ' 

where ao is the scale factor determined under the initial conditions. At £ >> A, 
a remains constant with increasing t. The other case that £ << A is likely not 
accessible, no matter which values will be assigned to a. Otherwise a slightly 
increases with increasing the comoving time t. 



3 Results and Discussion 

In Fig. [TJ H(t) and a(t) are depicted in dependence on the comoving time t. 
We compare our results, Eq. |TB] & UZl with H{t) and a(t) resulting in when the 
background matter is assumed to be filled with an ideal and non-viscous fluid 
and utilizing equations of state of non-interacting ideal gas 

H(T) = l/2f (18) 
a(t) = \Tt (19) 

In the case of viscous fluid, it is obvious that H(t) = a/ a has an exponential 
decay, whereas in the non-viscous case H(t) is decreasing according to Eq. [TH1 
The latter is much slower than the former reflecting the nature of the exponential 
and linear dependencies. The other difference between the two cases is obvious 
at small t. We notice a divergence or singularity associated with the ideal non- 
viscous fluid, Eq. [T8l The viscous fluid results in finite H even at vanishing t, 
Eq.QS 

Also the scale factor a(t) reflects differences in both cases. a(t) in a Universe 
with an ideal and non-viscous background matter depends on t according to 
Eq. 1191 which simply implies that a(t) is directly proportional to t and a(t) 
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vanishes at t = which reflects the singularity of H(t) = a/ a. Assuming that 
the background matter is filled out with a viscous fluid results in different a(t)- 
behavior with increasing t. At t = 0, a(t) remains finite. Correspondingly, H (t) 
remains also finite. In general, the dependence on t is much more complicated 
than in Eq. [121 Here we have an Aa root of an exponential function. If £ >> A, 
a remains constant with increasing t. 

As given in section[2l to deal with viscous Early Universe, one has to solve the 
second kind Abel first order ordinary differential equation [26] . which is a non- 
trivial task. In doing this, we have to approximate g(z) in the Abel canonical 
form, Eq. 1141 This is the only variable, which we approximated numerically. 

In doing this, we assumed that the Universe is flat, k = 0, and the back- 
ground matter has a finite viscosity coefficient. The resulting Universe is obvi- 
ously characterized by a constant scale factor and a vanishing Hubble parameter 
at large t. At small t, both a and H remain finite even at vanishing t, i.e., no sin- 
gularity. The validity of our solutions depends on the validity of the equations of 
states, Eq. [5J which we deduced from lattice QCD simulations at temperatures 
larger than T c w 0.19 GeV. Below T c , when the Universe cooled down, not only 
the degrees of freedom suddenly increase [29] but also the equations of state 
turn to be the ones characterizing the hadronic matter. Such a phase transition 
to hadronic matter would characterize one end of the validity of our solutions. 
The other end is at very high temperatures, at which the strong coupling a s 
entirely vanishes. 
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Figure 1: Top panel: Hubble parameter of viscous cosmology depicted as func- 
tion of the comoving time t, Eq. [T()]at a = — 1.0 and compared with the non- 
viscous case, H — 1 /(2£). Bottom panel: the scale factor vs. t, Eq.[l7]at a = — 1 
compared with the ideal behavior, a(t) = ^/t. 
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